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ABSTRACT 

We describe conformal field theories, correlation functions of which satisfy 
equations of the two-dimensional fluid mechanics. Prediction for the energy spec- 
trum is given, E(k) ~ A; -25 / 7 . 
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In this letter I will discuss the turbulent solutions of the equations of the two- 
dimensional fluid dynamics. The term "turbulent solution" means here that we 
will be dealing with the correlation functions of velocities, vorticities, etc. instead 
of these quantities themselves. Therefore, the implicit assumption in this paper 
is that in the turbulent flow we have some well defined stationary probability 
distribution with respect to which these correlation functions are defined. 

It is well known that, in this context, the Navier-Stokes equations lead to the 
chain of relations, expressing N-point function in terms of N+l point functions. 
The standard approach to the theory of turbulence is to use some closure hypothesis 
in order to get concrete equations for the correlators (see [1] for a review). 

Here we take a completely different route. We will try to satisfy the above 
mentioned chain of relations exactly by borrowing sets of correlation functions, 
known in the conformal field theory [2]. As a result we obtain a certain number of 
exact solutions for the two-dimensional turbulence. 

The basic equations, to be solved, are very simple. We have vorticity uj as our 
basic fluctuating variable. Navier-Stokes equations have the form: 

w + e a/3 d a il)dpd 2 il) = vd 2 uj (1) 

(Here d a = d/dx a ; d 2 = -j^. ip is a stream function related to vorticity by 
uj = d 2 i\) and to velocity by 

v is viscosity.) Now the correlation functions of the type 

< uj(xi) ■ ■ -u(xn) > 
must be time independent, which gives: 
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< u(xi)u(x2) •••> + < u(xi)u(x2) •••> H = (2) 

It is understood in (2) that we should express to through if; via eq.(l). This 
obvious equation (when being written for the generating functionals) is called Hopf 
equation [1]. It expresses the fact that the probability distribution of our system 
commutes with the equations of motion. Since in turbulence we are dealing with 
the large Reynolds numbers, the viscous term in (1) can be neglected. It plays the 
role of the ultraviolet cutoff for the theory and will be discussed later. 

Let us now present the simplest theory, which satisfies inviscid Hopf equation 
(2). Take the minimal model (2, 5) of conformal field theory (see [2] for the 
definitions). In this case we have two primary fields — the unit operator / and the 
field with the complex dimension (-1/5, -1/5), which we will identify with ip. The 
two point function of ip in the conformal theory is: 

< ip(x)ip(0) > = -|f| 4/5 (3) 

We shall postulate that, in order to get physical correlation function from the 
conformal one, we have to consider momentum representation of (3): 

< > = < 4 > 

and to transform it back by cutting off the /c-integral at some infrared point, defined 
by the large scales, k m i n ~ 1/r. This would give: 

< ^(£)^(0) > = const[R^ 5 - |£| 4 / 5 ]. (5) 

We will have in mind similar rule for all cases, namely that the momentum space 
correlation functions 

< tf)(ki) ■ ■■tj}{k N ) > 

for \kj\ >> 1/R are R-independent, while in the coordinate space one should 
include contributions from the infrared corners. Thus defined, physical correlators 
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differ from the conformal ones by the terms, analytic in xi — xj. It should also be 
understood that conformal expressions are valid only if 

— * 

\k\ < k max rsj l/a 

The ultraviolet cutoff a determined by the viscosity. The presence of the UV 
cutoff means that we should use careful point-splitting procedure when defining 
the products of operators. In particular the product in the LHS of (1) should be 
understood as: 

e a/3 da^(x)d/3d 2 ^(x) = \ime af 3d a ^(x + a/2)dpd 2 ip(x - a/2) (6) 

(where Urn means directional averaging). In order to evaluate the RHS of (6) 
we have to use operator product expansion, which, for our case, can be written 
symbolically as 

M x W\ = [i] + M 

where [ip] means conformal class of tp, i.e. itself together with the operators 

L- n being Virasoro generators [2]. More explicitly: 

ip(x + a/2)ijj(x-a/2) = \a\ 4/5 (l + Cia 2 L_ 2 / + • • •) 

+ \a\ 2/5 (ij + C 2 aL_ x i) + a 2 (C 3 L_ 2 + C 4 L 2 _^ + • • •) 
Performing differentiation of the LHS, we find the leading term: 

e al3 d a ^d d 2 ilj = cons*|a| 2/5 (L_ 2 Li 1 - L_ 2 L 2 _ 1 )V> (8) 

a — > o 

Antisymmetry of the RHS under complex conjugation is, of course, just the con- 
sequence of the e-tensor at the LHS. If we rescale 

i) => \a\ 2 '^ (9) 

requiring fixed values of velocity at the scale a, the a-dependence of the RHS 
disappears. It is now quite easy to see, that actually RHS of (8) is equal to zero. 
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This happens, because in the (2, 5) model, the field ip is degenerate on the second 
level, satisfying: 

(L_2 - 5/2Lii)^ = (10) 
Then, because of antisymmetry we get: 

u = -e a pdaipdf3d 2 i) = 

and the Hopf equation is formally satisfied. The only thing to check is that this 
result is not spoiled by the analytic terms, coming from the infrared contributions, 
mentioned above. Indeed, we dealt with the strictly conformal theory, disregarding 
these terms. In order to check this let us look at the expression for the operator 
(8) inserted into some correlation function in the momentum space. We have: 

< • • • iK/tf) > = Jd 2 k[k,q}(k 2 -(q-k) 2 ) 

x<^(k)^(q-k)^(h)--^(f N )>. 

The dangerous region is k — > 0. By using operator products once again it is easy 
to see that: 

<mi>{z-mh)-> k Zo]k-^ (12) 

and hence there is no infrared divergence in general. 

Let us estimate the region in which our solution works. First of all we compare 
kinematical time at the scale r: 

1 - ■ w(r > ~ ^(r)/r 2 (13) 



t(r) 

with the relaxation time: 



J- ~ M (14) 

r(r) ij;(r) 1 } 
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Using the fact that 



<Kr)~^(a)(-) 2 / 5 

(X 



and eq. (1) we obtain 



1 if}(a)a 2 



•1/r 2 



(15) 



r(r) r 



On the other hand, our neglect of the viscosity is justified if 



1 



> vjr 2 ; 



> 1 



(16) 



t(r) 



And hence 



a ~ R ■ (Reynolds number) 5 / 2 



(17) 



It should be realized that the solution we discussed is rather exotic. It breaks 
parity, since the stream functions are pseudoscalars and hence can exist only in 
the systems with net rotations or magnetic fields. Another unusual feature, which 
we will discuss now is the fact that as v — > the fluxes of energy and enstrophy 
vanish. This regime requires a subtle balance between injected rotation, enstrophy 
and energy. Still it may be interesting as a first case in which the Hopf equations 
are exactly solved (apart from the trivial Gibbs distributions). 

Viscosity and stirring forces should play the role of the boundary condistions 
for our solutions. At present we do not know any rigorous way to implement them, 
and will have to make some conjectures, which lead to more generic and physical 
solutions than the above. 

The standard picture of turbulence involves Kolmogorov's idea of constant 
fluxes. According to Kraichnan [3], in two dimensions the relevant flux is the flux 
of enstrophy. Perhaps the condition of constant flux is necessary for the matching 
of the perfect fluid solution with the viscous region. Let us try to formulate this 
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condition in our language. The enstrophy is a conserving quantity for the perfect 
fluid given by: 



H 



J co 2 (x)d 2 x (18) 



In order to study this quantity in the turbulent regime it is convenient to 
introduce the point-split definition (6) (which actually has been exploited already 
by Kolmogorov in 1941) and to use the identity 

d d 

— < u(x+r/2) > u(x-r/2) >= - — < u(x+r I2)u(x-r /2)(v a (x+r /2)-v a (x-r /2)) > 
dt or a 

(19) 

The constant r-independent enstrophy production will be reached if: 

< u(x + r/2)u(x - r/2)(v a (x + r/2) - v a (x - r/2)) >= Br a (20) 

This is a standard Kolmogorov-like relation. It has striking analogy with the axial 
(and other) anomalies in quantum field theory. Indeed, if one considers massless 
quantum electrodynamics, and define the axial current as 

J^{x) = lim < i)(x + r/2)7 M 75^(a; - r/2) > (21) 

the Dirac equation readily gives: 

df, J^{x) = lim(A M (a; + r/2) - A^x - r/2)) < ^j(x + r/2) lfll5 ^(x - r/2) > 
= const J e^Xpd^AvdxApd^x ^ 

(22) 

Physically, chiral charge is produced at the cut-off momenta and transferred through 
the momentum space into the physical region. In Kolmogorov's case the same hap- 
pens with enstrophy or with energy. 

While eq. (20) is very nice, the next standard step looks quite appalling to a 
field theorist. Namely, one concludes from (20) that dimensionality in this formula 
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can be counted as a simple sum of dimensions of 10 ~ v/r and v thus giving 
v ~ r. The problem here is that in the field theory dimensions are not additive 
and, when we take an operator product, we get at the right-hand side some defect 
of dimensions. 

We will attempt here to formulate constant flux condition field-theoretically. 
Let us consider some conformal field theory with the fusion rule: 

M x M = [0] + • • • (23) 

As was argued above we have 

Co ~ (L_ 2 £ 2 i - £-2£-i)0 (24) 
Now, if we rewrite eq. (19) in the form: 

< Co(x + r)uj{x) > = r° 
we obtain (recalling that uo = d 2 ijj): 

(A + 2) + (A^ + 1) = (25) 

A^ + A^ = -3 

If dimensions had been additive A^ = 2A^ (which is not the case of course) we 
would have obtained Kolmogorov-Kraichnan result Aw, = —1. We do not insist 
now on the especially slow relaxation, by requiring, as we did in our first example, 
that 4> is level two degenerate. Instead, we need the condition 

A > 2A^ (26) 

which simply means that RHS of (23) vanishes as a — > 0. In this case we get 
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inequality 



A^ < -1 



(27) 



which means that the energy spectrum, given by 



E(k) ~ k 



is steeper than in Kraichan-Kolmogorov approximation 



E(k) ~ AT 3 . 



Of course, it should be remembered when implementing (24) that non-zero con- 
tribution there comes from the regular (infrared) term, while the conformal term 
is zero, due to orthogonality between the primary operators. This is just as well, 
because otherwise logarithmic terms could be expected. Also, the conformal part 
of the theory is invariant under time reflection, while the correlator (24) breaks it 
due to "infrared leakage" of time asymmetry. 

I do not have complete classification of conformal theories, satisfying above 
requirement. There exist an appealing example, however. Let us consider minimal 
models of the type (2, 2 N+l). From our general classification [2] it follows that 
in this case we have a set of N primary fields, 



with dimensions: 



-A 



s 



(2N-s)(s-l) 
2(2iV+ 1) 



(29) 



and the algebra: 



[lfas-1 + [lp2s-3] H 



(30) 



for 2s — 1 < N, and with replacement s =>- 2N + 1 — s otherwise. 
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Constant flux condition (25) takes the form: 



A s + A 2s _i = -3 



(31) 



Solution to this Diofantian equation is s = 4 and iV = 10 for which case 



= A 4 = -8/7; A = A 7 
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y 



(32) 



The energy spectrum for this case is 



E(k) ~ l/k 3+4 / 7 



(33) 



The experimental value of this exponent seems to lie between 3 and 4. Also, this 
solution preserves parity, because there is no -^-operator in the product ip x ip, and 
more generally we have symmetry under parity reflection: 



There could be other solutions of these requirements. The major difficulty of 
our approach is the somewhat heuristic treatment of the analytic contributions. 
We need to supplement conformal field theory with more definite rules for their 
contributions. Before that, our analyses of the basic equations is not complete for 
large enough negative A (our first example with A = —1/5 seems to be immune 
to this problem). Probably, one possible approach is to perturb conformal theory 
away from criticality and create a small mass gap which will provide IR-cutoff. 
This is left for future work. 

The whole picture looks as following. There are many formal solutions of the 
inviscid Hopf equations. Some of these solutions (for which <fi is level two degen- 
erate) are more stable than the others. However, only one or very few solutions 
satisfy right boundary conditions (in the momentum space) needed to match them 
with the viscous region. All that is very much like the situation for the laminar 
flows, where one needs boundary layer consideration in order to determine the 
correct inviscid solution. 



V, a -(-l) a+1 ^ 



(34) 
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Another view on the same problem is the use of Lagrangian coordinates for the 
fluid. In this formulation there is some similarity to string theory which might be 
useful. Namely, it is easy to see that if we describe the fluid by the fluid coordinate 



= 1,2 



(35) 



VA(x(£,t),t) 



There exist a convenient action for the incompressible fluid (in any dimensions): 



where p is the Lagrange multiplier. 

This action has invariance under volume preserving diffeomorphisms, and should 
be treated in a stringy way, by introducing ghosts, auxilliary metric etc. Perhaps 
in the steady state (for D=2) we have some effective string theory describing cor- 
relations of x — s or, better to say, of the "vertex operators" : 



Finally, there are many obvious generalizations of the above approach, like includ- 
ing compressiblity, magnetic fields, passive scalars etc. 

I am grateful to E. Siggia, V. Yakhot and A. Zamolodchikov for useful discus- 
sions. I am also pleased to thank Kurt Gottfried and Andre LeClair for their kind 
hospitality at Cornell University, and for providing nice weather, which helped to 
complete this work. 




dtd D Za?(Z,t) + / p(#)[det(_±) - l]d\dt 



a 



(36) 
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